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NEWTONIAN GRAVITY OR GRAVITY ANOMALIES? 


Stacey and coworkers (Stacey and Tuck, 1981; Stacey et al. , 1981; 

Stacey, 1983; Holding and Tuck, 1984; Holding et al., 1986), and Hsui (1987) 
all present tentative evidence that gravity is non-Newtonian. Their geophy- 
sically determined G, the universal constant of gravitation, is on the order 
of 1 per cent larger than the laboratory value of 6.672 ± 0.004 x 10 -11 m 3 
kg -1 s“ 2 . Their values are derived from measurements of the gravitational 
acceleration g conducted in mine shaf ts, boreholes, and ocean depths, and a 
(Newtonian) theory of how g varies inside the Earth. 

Briefly, the technique is to assume a rotationally flattened Earth and 
to remove ellipsoidal shells above the depth at which g is measured. Since 
the gravitational attraction is zero inside a homogeneous ellipsoidal shell, 
only the ellipsoid below the gravimeter contributes to g. By correcting for 
the mass of the shells, one knows how much mass remains in the ellipsoid, 
and therefore its pull on the gravimeter. Departures from the expected 
Newtonian values at varying depths inside the Earth provide the evidence for 
non-Newtonian gravity. 

However, it is the local density that appears in the equations of the 
theory, and densities in the local environment are carefully measured. But 
the local density is not usually characteristic of the average density of 
the shell which must be used in computing its mass. For example, the local 
density of a borehole on a continent is about 2800 kg m -3 , while at a depth 
just below sea level the average density will be approximately (0.7) 1000 + 
(0.3) 2800 - 1540 kg m -3 , reflecting the various proportions of water and 
continental rock. This example illustrates that the average density can be 
much different from the local density. 

We wish to contnent on the resolution of this apparent contradiction. 

More importantly, we will demonstrate how large-scale density anomalies in 
the Earth (scale length > 500 km), if not corrected for, may masquerade as 
non-Newtonian gravity, and provide spurious evidence for a "fifth force" . 

We do this by showing that simplified Earth model containing only long-wave- 
length gravity anomalies can affect the measurement of G at the level 
reported for non-Newtonian behavior. 

In illustrating our remarks we assume the following: the Earth is non- 
rotating and spherical in shape with radius R and mass M. For simplicity in 
our expressions we shall not consider the centrifugal acceleration and 
flattening due to rotation because they can be analytically corrected for 
(Stacey et al., 1981; Dahlen, 1982). We further assume an Earth model whose 
density distribution p(r) has zero multipoles beyond a certain degree L but 
is otherwise identical to the Earth, so that all the sums appearing later 
are finite. In terms of spherical harmonics, p( r) can thus be written as 
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and Pj- is the normalized associated Legendre function of degree 1 and order 
m, ana r = (r,Q) where 0 is an abbreviation for both the co-latitude 6 and 
the longitude X. The spherical harmonics are orthogonal and we adopt the 
normalization commonly used in geophysics: 
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over the unit sphere. Apart from the rotational flattening, only the 1=0 
(monopole) terms have been considered previously. We shall consider the 
0 < 1 £ L terms as additional density anomalies . in the following we fix 
the origin at the Earth's center of mass so that 1=1 terms vanish. For 
1=2,4, rn=0 terms we consider only the departure of the Earth from the state 
of a hydrostatic flattening due to rotation (e.g. , Jeffreys, 1976). Also, 
ignoring any polar motion (of the rotational axis) which only produces a 
centrifugal acceleration on the order of 10 -7 m s~ 2 (or 10 /xgal) , we set 
1=2, m=l terms to zero. 

The gravitational potential at a point inside the Earth (r<R), satis- 
fying Poisson's equation V 2 U(r) = -47rGp(r) , is 
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Writing the first integral in equation (4) as J Q - J and considering the 
radial component by differentiation leads to ' 
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where the slight difference in the local vertical and radial directions _ is 
ignored. The C-u-i are the normalized multipoles of the density p( r), given 
by 
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where V is the volume of the (spherical) Earth. They arise as coefficients, 
often called the Stokes coefficients, of the external (r>R) gravitational 
field: 
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and have been routinely determined from geodetic measurements (e.g. , I^rch 
et al., 1981; Rapp and Cruz, 1986). 

We now differentiate equation (5) to compare with previous studies which 
measured the radial gradient of g(r) : 
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and gp = GM/R^ = 9.8209 m s“2 is an accurately determined Earth parameter. 
Note that the effect of the density anomalies inside the shell (of radius r) 
is absorbed in the C-^ coefficients (in so that only the density of the 
shell has to be explicitly considered (in P) . In finding Ag between two 
depths, as do Stacey et al. (1981), one integrates equation (8). This is 
why an integral involving the local density p( r) appears, and not the 
average shell density. 

Let us now examine^the influence of long-wavelength terms previously not 
considered. We write P(r) = Pp(r) + P(r) and Q(r) = Q 0 (r) + Q(r) where 
subscript 0 indicates the 1=0 (monopole) terms, and no subscript indicates 
the sum of all higher-degree terms (2 £ 1 £ L). From equation (8), the 
geophysical solution of G can then be written as 
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Previous studies only considered the rotation, the flattening (embodied in 
an 1=2, m=0 term), and the P Q (r) and Qn(r) terms in their solution of G. 

Our equation (11) constitutes a generalization of their formula to include 
P(r) and Q(r) which arise from global density anomalies, and reduces to it 
when P(r) and Q(r) are ignored. 

We now examine the correction terms P(r) and Q(r) . We will estimate 
their importance in equation (11) for the long-wavelength anomalies in the 
Earth, for which we have chosen L=36, corresponding to a scale length of 
>500 km. 

We consider P(r) first. If we assume the continents to have constant 
densities and vertical boundaries for the depths considered, then Pirni( r ) = 
constant and can be moved outside the integral in equation (10). Iflz << 
R, where z = R - rjls the depth, then this integral becomes approximately 
2(2l+l)z, so that P(r) - 2(z/R) p( r). For typical depths (e.g. , in a 
borehole) of 1 km, P(r) - 0.0003 p( r) and can be neglected to this order in 
equation (11). 

Under this approximation, the effect of the gravity anomalies on the 
solution of G, as a function of the geographical location, is 


AG(r) 2 g 0 Q(r) / 47rRp(r). 


( 12 ) 
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This is the error one makes in the solution of G if the gravity anomalies 
are ignored. Note that only the multipoles of the density distribution (the 
Ct_£ coefficients) are involved — the complete knowledge of the density 
itself is not necessary. 

We next evaluate Q(r). For the long-wavelength gravity anomalies we are 
considering (1 S 36), and depths z « R, the factor (R/r)^ 3 in equation (9) 
can be taken to be 1 with an error of ~(l+3)(z/R). Thus, for depths of ~1 
km, equation (9) may be evaluated at the Earth's surface with only about a 
0.6% error in Q(r) (not to be confused with the error in the G solution). 

For the Cy coefficients we use those of the GEM-10B gravitational field 
(Lerch et at., 1981). 

Using the resulting Q(r) in equation (12) gives <AG>/G £ 0.43%, where 
<AG> is the root-mean-square of AG(r) over the Earth's surface, with AG(r) 
being computed at the center of 5° x 5° squares. . The gravity anomaly- 
induced AG(r) can be positive or negative depending on where the measurement 
is made. Extreme values of AG(r)/G reach as high as ±2.5% for the GEM-10B 
field. These percentages are all in the range of the values reported, and 
demonstrate that long-wavelength anomalies can significantly affect the 
geophysical measurement of G. 

Stacey et al. (1981) have briefly discussed the free-air gravity 
gradient at the Australian mine site. They found a discrepancy of 8 X 10 _ ° 
sec -2 between the measured value and the computed value based on the mono- 
pole Earth model. Although as great as 10% of the gravity gradient measured 
in the mine shaft, it was absorbed into the uncertainty in their quoted 
final value for G. Here we wish to point out the following. The free-air 
gradient has the same form as our g 0 Q(r)/R (but evaluated outside the 
Earth), and hence the same physical source. The large discrepancy is thus 
an indication of the importance of this source, which we have identified to 
be the density anomalies — not just local but long-wavelength as well. 

This cannot be treated as an uncertainty but rather a systematic correction 
that must be made before a meaningful G value can be deduced. Note also 
that, at the Australian mine site, the sign of the free-air gradient discre- 
pancy is consistent with one that can cause a positive AG measurement as 
reported (Stacey et al., 1981). 

The modest aim of the present paper has been to indicate that even long- 
wavelength gravity anomalies may seriously affect gravity measurements, 
leading to an incorrect value of G. However, we make no " long-wavelength . 
corrections" for G for any particular location for the real Earth. This is 
because of the Kaula's rule-of-thumb (Kaula, 1967). which states that the 
magnitude of the C-^ coefficients decrease as 10" 5 l -2 . This empirical law 
has been shown by Rapp and Cruz (1986) to roughly hold to at least L=180. _ 
The sum in equation (9) as a result grows rapidly with 1 far past the limit 
of the GEM-10B field (L=36). Except for certain locations where the effects 
for IS 36 and 1>36 happen to have comparable magnitudes and opposite signs, 
the inclusion of higher degree terms for the real Earth will, in general, 
greatly augment the values based on GEM-10B. 
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